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Abstract. As part of a reanalysis of galactic Asymptotic Giant Branch (AGB) stars at infrared (IR) wavelengths, 
we discuss a sample (357) of carbon stars for which mass loss rates, near-IR photometry and distance estimates 
exist. For 252 sources we collected mid-IR fluxes from the MSX (6C) and the ISO-SWS catalogues. Most stars have 
spectral energy distributions up to 21 /tm, and some (1/3) up to 45 /jm. This wide wavelength coverage allows us 
to obtain reliable bolometric magnitudes. The properties of our sample are discussed with emphasis on ~ 70 stars 
with astrometric distances. We show that mid-IR fluxes are crucial to estimate the magnitude of stars with dusty 
envelopes. We construct HR diagrams and show that the luminosities agree fairly well with model predictions 
based on the Schwarzschild's criterion, contrary to what is widely argued in the literature. A problem with the 
brightness of C stars does not appear to exist. From the relative number of Mira and Semiregular C- variables, we 
argue that the switch between these classes is unlikely to be connected to thermal pulses. The relevance of the two 
populations varies with the evolution, with Miras dominating the final stages. We also analyze mass loss rates, 
which increase for increasing luminosity, but with a spread that probably results from a dependence on a number 
of parameters (like e.g. different stellar masses and different mechanisms powering stellar winds). Instead, mass 
loss rates are well monitored by IR colours, especially if extended to 20 /jm and beyond, where AGB envelopes 
behave like black bodies. From these colours the evolutionary status of various classes of C stars is discussed. 

Key words. Stars: mass-loss - Stars: AGB and post-AGB - Stars: carbon - Infrared: stars 



1. Introduction 

Winds from AGB stars replenish the Interstellar Medium 
with a large portion of the matter returned from stellar 
evolution (70% according to ISedlmavdl 19941) through the 



formation of cool envelopes JWinters etall I2002L 



where dust grains condense ijSchirrmacher et al.l 



2003) 



2003: 



We do not have yet a quantitative description of 
AGB winds, though attempts at m o delling them are 
a few decades old llSalpeti 

Il985l iGail fc SedlmavJ ll987| T More recently, hydro- 
dynamical and phenomenological studies of pulsating 
stellar atmospheres and of the associ ated mass loss 



ICarciofi et alTl2004h . Such grains carry the signature of 
the nucleosynthesis episodes occurring in the AGB phases 
ijBusso et all llQOflt IWasserburg et al. 1 1200.^ : their pres- 
ence has been recognized in meteorites and brings direct 
informatio n on circumstellar processes (e.g. IZinne 

2001J). As AGB stars radiate most of their flux at 
long wavelengths, large surveys of infrared (IR) observa- 
tions play a fundamental role in studying their lumino sity 
and their winds (see e.g. lHabmdll99f)[ lEpchteinl ll999i. 
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Fig. 1. SWS spectra for an irregular variable (TX Psc), a 
Mira variable (IZ Peg) and a post-AGB star (V353 Aur). 
The dominant role of IR emission longward of 20 /im for 
Miras and post-AGB sources is a general property of our 
sample. 



proved knowledge o f stel lar distances (e.g. IVan Eck et all 
FflflilKnann e t al.l l2003t Iplrgea.t fc ChevaHieHteOOfft . 

A real step forward would be to derive realistic for- 
mulae, linking the efficiency of stellar winds to the lumi- 
nosities, colours and chemical properties of AGB stars, to 
be adopted as inputs for stellar models, thus avoiding free 
parameterizations. Thi s kind of studies has b ecome quan- 
titative only re cently llvan Lo on et al. l 12005ft . while pre- 
vious attempts ijVassiliadis fc Woodlll993t iBlockerlllS 
suffered from large uncertainties l|Woodlll996j) . 

Reducing those previous uncertainties also requires 
a good knowledge of absolute magnitudes, through im- 
proved distances; so far, mass loss studies for galactic 
AGB stars often ad opted some average constant value for 
the luminosity (see lJural Il986t lJura k, Kleinmannlli989l 
iLe Bertre etallhoOlir Relevant except ons exist, e.g. in 



the thorough a nalysis of C-rich source s by Bergeat and 
his group (see iBergeat et al 1 l2002alb! . hereafter BKR1, 
BKR2), where howeve r the photometric data are rather 
heter ogeneous (see also lKnapik et al.lll999t iBergeat et alJ 

With the above difficulties in mind, we plan to perform 
a thorough reanalysis of AGB luminosities and mass loss: 
this is actually the main scientific scope of our project 
for putting an IR telescope (the International Robotic 
Antarctic Infrared Teles cope) in Antarcti ca, at the Italian- 
French base of Dome C ijTosti et al 120041) . It will allow ex- 
tensive monitoring of AGB sources in Magellanic Clouds 
and in the southern Milky Way up to 20 /im (and possibly 
beyond, given the unique characteristics of the Antarctic 
atmosphere) . This paper is in fact part of the preliminary 
works necessary to define the key projects for the tele- 
scope. 

We shall therefore consider, in a series of works, all 
types of AGB stars (M, MS-S, C giants), using existing 
catalogues of IR observations and compiling a homoge- 
neous list of luminosities, distances, mass loss rates, look- 
ing for correlations between them and the IR data. We 
start with C-rich AGB stars; next steps will analyze MS-S 
giants, where C- and s-element enrichment becomes ob- 
servable, and finally the more disperse family of M giants. 

In this note we collected a sample of 357 C-rich sources 
with near-infrared (IR) photometry and mass loss esti- 
mates; these last are revised by scaling previous obser- 
vations with updated distance determinations. For an- 
other sample of 252 C stars we have collected also space- 
borne mid-IR photometry in the A = 8 — 15/Lim inter- 
val and/or at A ~ 21/xm. The intersection of these two 
groups is made of 214 C-rich objects. For this sample, 
near IR data are provided by ground-ba sed surveys like 
2MASS llCntri et al .1120031) and DENIS l|Er>chtein et a.lJ 
1999l lFouaue et al.l2000|) : those for mid-IR are taken from 
the recent space-borne missions, ISO and MSX. The rele- 
vance of MSX colou rs for classifying A GB stars has been 
recently pointed out (jOrtiz et al. 12005^ , with an approach 
which is complementary to ours. IRAS measurements are 
not included in our analysis because of their lower quality 
and spatial resolution (hence larger contamination) com- 
pared to more modern experiments. A homogeneous com- 
pilation of IR colours for stars of known distance will allow 
us to precisely estimate the bolometric corrections (BCs) 
and to correlate mass loss rates with luminosity (and/or 
the colour excess). 

We started from previous work done by BKR1 & 
BKR2, who classified C stars in several discrete families: 
the more advanced subclasses, named 'Carbon Variables' 
or CVs, from CV3 on, correspond to bright thermally 
pulsing AGB stars. Bolometric magnitudes were esti- 
mated by thos e authors using op tical photom e try b y 
iBaumertl l|l972j) . near-IR data from iGezari et alJ l|l993h . 
and including IRAS fluxes at 12 and 25 /im. Careful 
upgrades of Hipparcos were included. In our work we 
shall consider newer and more homogeneous IR data 
from the ground and from space, with Spectral Energy 
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Distributions (SEDs) covering a wider wavelength baseline 
up to 45 /im. In this way we also aim at verify ing on wider 
statist i cs suggestions p r eviously adv a nced bvlBusso et alJ 
(|l996f) . iMareneo et all l|l997l Il999j) . iBusso et alJ l|200lj) . 
according to which mid-IR colours are good indicators of 
the mass loss efficiency, and permit a first classification of 
the chemical properties of the circumstellar envelopes. 

The choice of discussing first a sample of Carbon- 
rich AGB stars is motivated by the fact that optical 
C(N) stars (usually irregular or semi- regular pulsators) 
present a sufficient level of homogeneity to be consid- 
ered as a snapshot on the relatively long (1-3 xlO 6 yr) 
TP-AGB evolution. They correspond to the moment in 
which the C/O number ratio reac hes unity in stars of 
moderate mass (l o wer than 3 IcTaussen et al.||l987t 
I Busso et all Il999t lAbia et alJ boOll 120021: iKahane et all 
2000). When looking at IR properties, however, one has to 
remember that more massive and evolved sources become 
part of the sample, thus ma king the picture more com- 
plex l|Barnbaum et alJll99l|) . Mass loss rates have been 
shown to span a wide interval, from a few 10~ 8 M Q /yr, 
l|01ofsson et alll 993alhl:ISchoier fc Olofssor^OoTl) to 10" 4 
Mfp/yr for the most mass ive and/or evolved objects 
llOroenewegen et a,l.ll2002alhh. 

All stars outside the mass range in which C-rich at- 
mospheres are formed will reach the final stages while re- 
maining O-rich. This is due either to the inefficiency of 
dredge-up (for lower masses), or to the fact that mas- 
sive AGB stars (M > 5 M Q ) burn the new carbon as 
soon as it is dredged to the envelo pe, thanks to the hot 
botto m burning (hbb) phenomenon l)Karakas fc Lattanziol 
l2004l) . Observational e vidence of this occurrence is emerg- 
ing in various cont exts (Smith e^ajL^^^^5; van L oon et alJ 
I2001L 1200.4 Iwhitelock et al]l200fl ICioni et a,l]l200,4 See 
in particular Ivan Loon et alJ l|l999albft . hereafter VL1, 
VL2. This paper is organized as follows. In section 2 the 
C-star sample and the photometric measurements are pre- 
sented; in section 3 we discuss mass loss rates from radio 
studies, together with the required updates (e.g. on the 
distance) . The luminosities and colours of our sample stars 
are presented in sections 4 and 5, together with the conse- 
quences of our analysis for the AGB evolution and dredge 
up. Correlations found between luminosity, colours and 
mass loss rates are reviewed in sections 6. Some general 
conclusions are drawn in section 7. 

2. C-star photometry from ISO-SWS and MSX 
catalogues 

Ground-based photometric studies of C stars in the 8- 
14 fim window were often performed through 1/zm-wide 
filters, centered around 8.8, 11.7 and 12.5 /im. We decided 
therefore to compile colours including filters as close as 
possible to the above ones, to make comparisons easier. 

The largest inventories of IR observations longward 
of 2.2 inn, for AGB s t ars co me from space experiments. 
iLe Bertre etall l|200ll |2003|) examined the sample from 
the Japanese experiment IRTS. We shall consider in- 



stead data from the European ISO mission and from the 
American MSX satellite (using the Infrared Point Source 
Catalogue, version 6C). In the case of ISO, the best known, 
nearby AGB stars are always too bright to be observed 
by th e imaging camera (ISOCAM: see ICesarskv et alJ 
Il996l) . Luckily, the database of the Sho rt Wavelength 
Spectrometer (SWS: FValentiin et al.lll99fih contains spec- 
tra of many known AGBs, and offers also a good photo- 
metric calibration. We therefore performed a re-binning 
of AGB SWS spectra from the ISO archives at the Centre 
de Donnees Stellaires ( CDS) of Strasbourg. In doing this, 
the spectral energy distributions were convolved with the 
response curves of common filters centered at 8.8, 11.7 
and 12.5 /xm, av ailable from the TIRCAM consortium 
l|Persi et alJll994j) . For the sake of simplicity, these filters 
will be indicated by: [8.8], [11.7], [12.5]. We also rebinned 
data at longer wavelength, convolving them with the fil- 
ters at 14.7 and 21.3 /im from MSX. Sources for which 
ISO fluxes in different epochs vary significantly were com- 
bined to get average spectra: a few sources showed the 
remarkable property of a global (bolometric) variation up 
to a factor of two. Such phenomena have been observe d 
previously l|Le Bertrelll99a Il993t Ivan Loon et"aflll998h: 
for int erpretations the reader is referred to lvan Loon et alJ 

dun. 

A number of our sources have complete SWS spectra 
up to 45 fim: those among them that also have near- 
IR photometry and estimates for mass loss and distance 
(49 in total of the about 70 astrometric sources) repre- 
sent a 'best selection', from which we can expect to derive 
all photometric properties with high precision, including 
the bolometric magnitudes (these last will simply come 
out from integrating the flux from 1 to 45 /im, as for the 
whole sample optical fluxes are very small compared to 
the IR ones). Three examples of such complete spectra, 
representative of stars with small, moderate and large IR 
excess, are presented in Fig.^ They show one of the main 
points we want to underline in this work: deriving bolo- 
metric magnitudes from too short wavelengths can be very 
misleading, especially for very re d sources. As examples , 
using BCs derived from J-K, as in lCosta k. Frogell l)l996jh 
the sources in Fig. ^ would be assigned absolute magni- 
tudes of: -4.8 (TX Psc); > 1.0 (IZ Peg); > 1.0 (V353 
Aur). The last two numbers derive from assuming gener- 
ous upper limits of 0.05 Jy for the very low J fluxes of 
these red sources. When BCs are instead computed start- 
ing from IR phot ometry extending from the K-band to 
10-20 /xm, as in 



ILe Bertre et al 



(2001), more reason- 
able magnitudes are obtained: Mboi — —5.2, —4.8, —5.9, 
respectively. When, finally we apply our BCs (see next 
sections) derived from photometry up to 45 /im, we get 
Mboi — —5.5, —4.8, —6.2 in the three cases. For very red 
sources (about 25% of Miras and most post- AGBs) even 
our Spectral Energy Distributions (SEDs) up to 45 /im 
are insufficient to gather the whole flux. For them, the 
luminosities we estimate will be lower limits. 

Concerning the MSX sample, we adopted the fluxes at 
14.7 and 21.3 /im as provided by the catalogue, while in 
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the 8-14 /im range, fluxes in the [8.8] and [12.5] filters 
were inferred from those in the MSX bands 'A' and 'C. 
In order to do this we selected a group of sources in com- 
mon between ISO-SWS and MSX. For them, we derived 
correlations between the two sets of filters, obtaining the 
following relations: 

logFi2. 5 = 0.80176 log.F c + 0.3349 



logFg.g = 0.9377 log Fa + 0.1266 

with regression coefficients (i? 2 ) of 0.95, 0.94 respectively. 
The relations apply to all sources examined: colour effects 
on different types of stars are negligible. They were used 
for expressing the MSX fluxes in the set of filters adopted. 
This can be done without loss of accuracy, as the uncer- 
tainty of IR colours is at least ± 0.1 m . No flux at 11.7 /im 
can obviously be derived for MSX stars. 

The resulting list of ISO and MSX C-rich AGB sources, 
with their inferred IR fluxes, is presented in Tables ^ to 
0| Here variability types are taken from the Combine d 
General Catalogue of Variable Stars l|Samus et alJl20o3l . 
In the tables we have included near-IR data from 2MASS, 
together with general information available on the vari- 
ability type and on parameters relevant for estimating 
mass loss rates (see next section). One problem affecting 
our sample is that data from catalogues are single-epoch 
measurements; though the amplitude of the light curves 
decreases sharply when moving from optical to near-IR 
wavelengths, we expect that J, H, K data can be still 
affected by uncertainties of several tenths of a magni- 
tude. This problem affects only the bluest semi-regular 
and irregular variables (at the lowest luminosities and 
mass loss rates considered). Indeed, despite the fact that 
redder sources also show larger amplitude variations in 
near-IR filters, their flux is usually dominated by mid-IR 
wavelengths, where the variability vanishes because the 
photosphere is not seen. The very large near-to-mid IR 
colours of several Miras and of all post- AGB sources make 
any variability at short wavelengths irrelevant (exceptions 
might occur for the few sources showing factor-of-two vari- 
ations in the bolomctric luminosities). Also, the photomet- 
ric data in the tables do not contain corrections for inter- 
stellar extinction, for two reasons. On one side extinction 
becomes small at the long wavelengths considered here; on 
the other hand, our stars are distributed in all directions of 
the sky, so that we should adopt average extinction laws, 
but these are rather uncertain, and it is not clear that 
(at IR wavelengths) average corrections are substantially 
better than nothing. As our stars are occasionally very 
distant (above 1 kpc) one has to estimate the effects of 
neglecting extinction. Let's consider representative cases 
at 200 and 2000 pc distance. For the average interstel- 
lar m edium one has A v (200yc) = 0.15 m ; A y (2000pc) = 
2.2 m llCillett et a,lJll975t iKnudelll 979ft. The w avelength 
dependence can be taken from iDraina l|l989l L as being 
A\ ~ A -175 . At mid-IR wavelengths (10 /im and beyond) 
this implies negligible corrections (0.015™ at maximum), 



well within the photometric uncertainly, even for the most 
distant sources. In the K-band (2.2 fim) extinction is ir- 
relevant within lkpc, but may be appreciable for the most 
distant stars (0.2™ 1 at 2 kpc). We notice that the unknown 
variability induces a scatter whose amplitude is similar to 
the uncertainty related to extinction. Neglecting this last, 
however, introduces systematic effects on the most distant 
stars, so that their near-IR luminosity will be underesti- 
mated. Hence at near-IR wavelengths, catalogue data for 
distant (> 1 kpc) LPVs taken at face values are typi- 
cally uncertain by a few tenths of a magnitude (from both 
statistical and systematic errors) . The effects on the bolo- 
mctric magnitude distribution in our sample are however 
small, due to the large IR excesses of most distant sources 
that make near-IR data relatively unimportant for them. 

When using the data of Tables ^ to 0] for determin- 
ing standard colours, we applied st andard photometric 
calibrat ions, adoptin g suggestions bv lBessell et al.l 1 1998?) 
and by lOlassl l|l999l) and interpolating them with black 
bodies for obtaining the fluxes in the selected l-/im-wide 
filters. According to this procedure, the fluxes (Jy) of a 
zero-magnitude star in the filters used are: 52.23 ([8.8]), 
29.55 ([11.7]), 25.88 ([12.5]), 20.25 ([14.7]), 8.91 ([21.3]) 
For near-IR, 2MASS calibrations are given in Cohen et al. 
(200§. 



3. Estimates for mass loss and related problems 

The work by l.lural l|l98fil) and l.lura fc Kleinmannl l|l989h 

set the stage for subsequent studies on AGB winds. These 
authors used mass loss rates estimated across the previous 
ten years, and mid-IR photometric observations from the 
IRAS satellite. Unfortunately, the scarcity of data then 
available led them to adopt average hypotheses on the 
luminosities, assuming a constant value of 10 Lq for all 
types of AGB stars. This fact (and the modest precision of 
IRAS photometry) prevents us from making too system- 
atic comparisons of those works with the present analysis, 
where the luminosities are estimated thanks to measure- 
ments of the distance, sometimes of astrometric quality. 
We shall however comment later on mass loss rates for 
the sources we have in common with the quoted studies. 

The inventory of more recent analyses on mass loss 
in AGB s tars includes works d e dicated to Carbon- 
rich stars ( Olofss on et al.l Il993albr iKastner et al. Ill993t 
Schoier fe Olofsson! l2001i IGroeneweeen et al.l l2002albt 



Schoier et all l2002t). to Oxygen-rich M-MS-S stars 
llSahai fc Liechtil Il995t IGroeneweeen fc de Jond Il998t 
lOlofsson et all20f)2ft . and more general approaches includ- 
ing a ll classes of AGB objects teastnerlfl992l lLoup et all 
| 1993[ iLe Bertre et"aH l200lL l2003t IWinters et al.l I2001 
2003). From the theoretical point of view, the very lim- 
ited evidence of coronal X-ray emission lljura fc Helflandl 
119841 ISahai et a,ljl2003t IKastner fc Sokerll2003h demon- 
strated that mass loss mechanisms were different than in 
the sun. Radiation pressure on dust grai ns was suggeste d 
to be at the origin of red giant winds bv lSakoeterl l)l974j) . 
In this approach, if the momentum of the radiation field 
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can be completely transferred to dust and gas particles in 
the envelope, then the mass loss rate can be expressed as: 

M ~ 2 x KT 8 T d (X*/L Q ) x (v/kms- 1 )- 1 M /yr 

where Td is the average optical depth of dust. The gas 
phase of the envelope might be dragged by dust, but a 
problem is how to transfer momentum to dust close tc 
the photosphere, as above 1500-2000 K no solid com- 
pound made with light elements (C to Si) would avoid 
evaporation. Hence, stellar puls ation in different modes 
l|Wood et alJll999t IWoodl l2000|) . pushing the radius tc 
large distances and the temperature to low values, and in- 
trod ucing shock waves should also be very important (see 
e.g. [Fleischer et alJ Il992t Il995t IWinters et all Il994albl . 
ll997ll2000aM . 

All approaches suffer from many limitations, so that re- 
cent work was concentrated on semi-empi rical approaches 
derived from observations. As an example . iLe Bertre et all 
l|200lL l2003h used near-lR fluxes (in K, L'), then adopt- 
ing simplifying assumptions like e.g a constant luminosity 
of 8000 Lq for all the sources studied. Being based on 
a wavelength interval where the AGB photosphere dom- 
inates the flux, these estimates are also affected by the 
intrinsic uncertainty of the variability. 

Other mass loss estimates were derived from obser- 
vations at millimeter or radio wavelengths, where the 
variable star is not seen. In most cases they adopt 
procedures early suggested b v iKnap p fc Mo rris] JT 985') 
for CO lines. For example. l Olofsson et alJ (J1993a b) 
and iGroenewegen et al.1 l)2002albl) assume emission from 
an optically thick environment, with a fixed size for 
the CO envelope, Rco- Due to this, they can some- 
times underestimate the wind efficiency by factors 3-4 
^Scnoier fc OlofssoulEiol. 

A m ore gener al case w as considered b v lLoup et al.1 
(119931). iKastnerl lll992|) . iKastner et all (|l993j) and 
IWinters et al.1 l|2003j) . The results are given in terms of 
the distance d (in kpc), of the expa nsion velocity V K ( in 
km/sec), and of the flux /. We adopt "Loud et al.1 l|l993h 's 
formulation whenever possible, re-scaling their mass loss 
rates with more rece nt estimates for th e distance and for 
V e . When data from lLoup et alJ lll993l) are n ot available, 
we use instead the lOlofsson et al.1 <|1993albl) procedure, 
again updating the parameters when possible. 

Tables H to 0] report the relevant data for those stars 
for which IR photometry was derived in the previous sec- 
tion, together with the updated estimates for mass loss 
rates. It is interesting to note (see Fi g:. that fo r the 
few sources we have in common with I Jural (|l986h and 
I.Tura & Kleinmaiml l|l989|) and despite the mentioned un- 
certainties of these older analyses, there is a fair correla- 
tion between the mass loss rates given by the above au- 
thors and those derived by us (with a regression coefficient 
R 2 = 0.88). In most cases the sources are Mira variables. 

We shall be guided by the data of the AGB stars 
(Table ^) with the best distance estimates. They in- 
clude objects with astrometric measurements (usually 
from Hipparcos), for which we adopt the recommended 
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— y = - 0.51 +0.85 x 



dM/dt (J86.JK89) 

Fig. 2. The existin g corr e lation betw een mass loss esti- 
mates published by ljural l|l986|) and lJura fc Kleinmanr] 
jl9_8_9]) and our re-calibration of measurements from 
lLoup et ail 0*993) with updated distances. 

values bv lBergeat fc Chevallierl l|2005h . plus a small group 
of post- AGB stars for which the distance was derived from 
dedicated studies of the circumstellar envelopes. Our men- 
tioned 'best selection' is actually made by those sources, 
among the objects in Tabled for which also the integral of 
the SED up to 45 /im is available. We give a lower weight 
to other methods of deriving distances (the corresponding 
points in the figures will be smaller). When nothing else 
is available, we qu ote in the tables distances estimated by 
lLoup et al.Nl993h by assuming a specific value for the lu- 
minosity. Obviously, these sources will be used only when 
a precise knowledge of the distance is irrelevant, e.g. in 
colour-colour diagrams. 

4. Energy distributions and photometric 
calibrations 

A few SEDs, chosen to represent the various shapes found 
in our sample, are shown in Fig. They illustrate how 
different the spectra for different types of sources are, 
displaying the increasing importance of mid-IR (10 and 
20 /xm) wavelengths when moving from irregular (or semi- 
regular) pulsators toward Miras and post- AGB stars. This 
has two important consequences as discussed below. 

First of all, it is clear from the figure that, when the 
observations are not extended up to at least 20 /im, most 
of the flux from the reddest objects is missed. Plots like 
those in Fig. ^ reveal that many AGB stars radiate the 
majority of their flux in a region between 20 and 45 /urn. 
This range is rarely considered in dealing with large sam- 
ples of sources. In the past, the flux in the 60 /im IRAS 
filter was often used as representative of far IR, and corre- 
lated with mass loss. However, IRAS 60 fim data may be 
partially contaminated by diffuse (cirrus) emission; more- 
over, the IRAS spatial resolution at those wavelengths is 
very poor. 
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Fig. 3. The Spectral energy distribution up to 21.3 /jm 
of a number of sources in our sample, divided according 
to their type, as constructed with the set of l-/^m-wide 
filters described in the text. 
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Fig. 4. Model bolometric magnitude of a 2 M & AGB star 
during the phase of thermal pulses on t he AGB. Model 
from t he FRANEC code, as computed bv lStraniero et all 
(1991. 



Secondly, when we consider the relative number of 
sources having SEDs with the various appearances of Fig. 
(with semiregulars and Miras occurring with nearly 
equal frequency) doubts arise on the common idea that 
the different variability types are alternatively encoun- 
tered along the AGB due to the variations o f luminosity 
in co rrespondence of thermal pulses (see e.g. ICioni et alJ 
2001). These doubts are best illustrated with reference 
to the behaviour of the stellar luminosity on the AGB 
(Fig. UJ. The usual interpretation might in principle be 
true for the first thermal pulses, when the star is oxygen- 
rich, as regimes of higher and lower luminosity are both 
encountered, due to the occurrence of post-flash dips. In 
early phases these dips cover nearly 30% of the time be- 
tween adjacent pulses (for a recent review on the models, 
see e.g. lStraniero et al.ll2005|) and hence imply that lower- 
luminosity variables (identified with semiregulars) should 
account for ~ 30% of the total number of O-rich AGB 
stars. However, for the final C-rich phases the duration of 
low-luminosity regimes becomes increasingly short with 
increasing pulse number. According to Fig. relative to 
a 2 Mq AGB star of solar metallicity, the stages where the 
luminosity is lower than average by at least 0.2 m account 
for about 10 percent of the time spent by the star as a 
C-rich giant: if C-rich semiregular variables were to be as- 
sociated with these low luminosity regimes, they should 
be pretty rare, contrary to the evidence we have. Our 
sample was selected from mass loss data and therefore is 
certainly biased toward the reddest and more mass-losing 
Miras, but despite this semiregular variables are of equal 
statistical weight. If the AGB stars really switch from one 



variability type to another (see e.g. lCioni et al.ll200l[) this 
is unlikely to be related to thermal pulses; rather, in atmo- 
spheres animated by complex pulsations, superposition of 
close-by frequencies can induce beats and hence amplitude 
modulations (with p eriods of centuries or millenia, see e.g. 
iMareneo et al.lE00l|) that may be independent from what 
happens in the internal layers. One has to notice that, in 
the BKR1 & BKR2 classification in CV classes, the com- 
plex behaviour of the C/O ratio shows an increase with 
class number and with decreasing T e ff, down to about 
2500K; Miras (classes CV5— CV6) occupy regions with 
higher C/O ratios than SRa,b stars (classes CV3+ and 
CV4). This fact, and the gradual SED changes we see (Fig. 
13), with smooth transitions from Semiregulars to Miras 
and then to post-AGBs, leads to the conclusion that, at 
least on average, most Miras are in an evolutionary stage 
subsequent to semiregulars. Oscillations between the two 
types should exist, but the final evolutionary stages should 
be more heavily populated by Mira-like variables. 

SEDs like those of Figure 01 were the base for esti- 
mating bolometric magnitudes. We used for this the fun- 
damental re lation of p hotometry, according to which one 
has (see e.g. lGlasdll999|) : 



Mi 



bol 



F u dv + C 



(1) 



where C = —18.98 when the flux is expressed in W/m 2 . 
As a consistency check we derived the zero-magnitude flux 
from the data of iBessell et al.l l|l998h . extended from the 
U to the Q band, finding C=— 19.01, a very good approx- 
imation to (1). The bolometric magnitude of AGB stars 
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Fig. 5. Panel a): Bolomctric corrections for the K magni- 
tude as a function of the K-[12.5] colour. They were de- 
rived for AGB C-stars with complete SEDs, from 2MASS 
and ISO-SWS, up to 45 fim. Panels b) and c): bolomet- 
ric corrections for the [8.8] magnitude, as a function of the 
mid-IR colour [8.8]-[21.3]. Of the two relations shown, the 
one in panel c) is more appropriate for extrapolations to 
very red sources, as the cubic spline of panel b) introduces 
unwanted biases outside the range shown. 

has often been derived from corrections applied to the K 
magnitude; these corrections are usually a function of the 
K-[1 2 j colo ur (using e.g. IRAS data, see iLe Bertre et all 
l200ll 12003^) . For the sake of comparison, we too present 
here our Mj D j estimates through bolometric corrections as 
a function of the (K-[12.5]) colour. If / is a measure of 
the flux in Jy (Tables [2 Ej): 

K = -2.5 x log U, K + 5-5x log d + 2.5 x log f v , Ko (2) 

where d is the distance in parsec and j v ,K is the zero- 
magnitude flux in K. 

Bolomctric corrections derived from equations (1) and 
(2) are shown in Fig. (panel a) as a function of the K- 
[12.5] colour, together with a fitting formula for extrapola- 
tions. Optically-selected sources (in general semi-regular 
variables) have small corrections: their flux is well esti- 
mated from traditional criteria at short wavelengths, as 
most flux is radiated in near-IR. Alternative expressions 
of the B.C., as a function of mid-IR colours alone, to be 
used for very red sources too faint in near-IR, are pre- 
sented in Fig. (panels b, c). Here semiregular variables 
are omitted, because of their small IR excess (see Fig.^J. 

Finally, we must notice that mid-IR emission from cir- 
cumstellar envelopes makes the stars of our sample ap- 
pear as black bodies at long wavelengths, apart from the 



Fig. 6. Examples of colour-colour diagrams in the IR for 
C-rich AGB and post AGB stars. 

known typical features of C-rich dust. This is shown in 
Fig. where calibrated colour-colour diagrams are shown 
for near- and mid-IR filters. In the top panel, adopt- 
ing the K— [21.3] colour as baseline, the data points of 
known AGBs and of sources (mainly from MSX) of un- 
known variability type are almost perfectly aligned on the 
blackbody sequence; they should be dominated by cool, 
evolved Mira stars with extended dusty envelopes; post- 
AGB sources instead display excess emission at 21.3 /im, 
which is a well-known proper ty of such C-rich objects 
(van Winckel & Revniersl UpOO) . The bottom panel then 
shows that, if only mid-IR is used, including the [11.7] fil- 
ter where SiC and PAH features are present, then normal 
AGB C stars show excesses in this filter, while post-AGB 
C-rich sources gradually move to an excess at 21.3 /xm, 
while becoming redder. In general, Miras are redder than 
semiregulars and Post-AGBs are the reddest sources of 
the sample (due to their cold dust). In the evolutionary 
hypothesis we have tentatively advanced before (at least 
on average), this graph would offer an immediate tool for 
classifying the evolutionary status of C-rich circumstellar 
envelopes. 

5. Bolometric magnitudes, colours and model 
comparisons 

The magnitudes derived by our analysis are used in Fig. 
where HR diagrams are presented, using an IR colour 
as abscissa. As stated, small symbols refer to stars with 
less-accurate distances: they are usually distant post- 
AGB objects, or AGBs of unknown variability type. The 
dashed area in Fig. [7| represents the zone covered by 
standard stellar models with the Schwarzschild's crite- 
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Fig. 7. The HR diagram of observed C stars in near IR, 
as compared to the area (dashed) covered by canonical 
stellar models (without hbb). The minimum limit for C 
star occurrence in the adopted models is indicated 



rion for convection, adopting the usu al mixing length the - 
ory as adapted to stellar physics bv lCox fc Giulil l|l968l) . 
The value of the so-called parameter a (ratio of the 
mixing length to the pressure scale height) is a = 2.1 
l|Straniero et alJ Il998h . All models considered are from 
the FRANEC evolutionary code; they refer to metal- 
licities from Z=Zq/3 to Z=Z a nd include models of 
relatively low ma s s (2 to 3 Mm) llStraniero et al.lll997t 
iGallino et al. I ll998t iBusso et al.lll 999). plus intermediate 
mass star models l|Dommguez et alJll999t IStraniero et alJ 
2000|. As representative cases, we plot explicitly the track 
of a 2M Q star with Z=Z Q , and that of a 3 Mq star 
with Z~Zq/3 (heavy lines). As compared to the evolu- 
tionary tracks, the data points (which are systematically 
displaced at redder colours) show the gradual formation 
of cool and dusty circumstellar environments, which recy- 
cle part of the stellar flux at longer wavelengths through 
the absorption and re-emission mechanisms of dust grains. 
For the bluest sources, the displacement in J-K with re- 
spect to theoretical predictions might also be due to mod- 
els underestimating molecular opacities of C-rich material 
l|Marigo et alJl2003l) . 

We underline that the models adopted do not as- 
sume any overshoot fo r convection (contrary to, e.g. 
iFrost fc Lattanziolll996|k for the techniques used t o com - 
pute the convective border see IStraniero et alJ l|2Q03(l . 
Overshooting from the convective envelope and in gen- 
eral an extension of dredge up larger than allowed by 
the Schwa rzschild's crite rion is even today a common as- 
sumption llHerwidEiol . This is usually claimed to be 
needed for letting the photospheres become C-rich at low 
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Fig. 8. Top panel: the histogram showing the fractional 
number of C stars per magnitude interval in our sam- 
ple (230 sources). Bottom panel: same as before, for 
the samples of LMC C stars in optical wavelengths 
ilCosta fc Frogellll99(l 887 sources) and in the infrared 
(VL2, 29 sources; note that the bin size here is the same 
as used by these authors for their data. 



luminosities, in order to fulf i ll con strain ts like those by 
iBlanco. Blanco fc McCarthvl l|l980h and ICosta fc Frogel 
(1993), derived by optical and near-IR data, suggesting 
average C-star magnitudes below —5. We notice that sim- 
ilar works in the optical bands became recently feasible 
on other galaxies, and the results were essentially the 
same, yielding C-star mag nitudes between —4.5 and —5 , 
again neglecting IR fluxes l|Battinelli fc Demersll20f)4albT : 
iDemers et aT~ll2004|) . 

Despite some contrary evidence, which has emerged 
fro m IR studies of Magell anic Clouds and of the Galaxy 
(see lvan Loon et al.ll|l998h . VL1, VL2, BKR2) some theo- 
rists assume, even today, that such faint bolometric mag- 
nitudes characterize the whole family of C-stars, and 
yield a disagreement wit h stellar models adop t ing the 
Schw arzschild's criterion jStancliffe et alJ l2f)05l iHerwiel 
2005). Our data do not support this belief: Fig. re- 
veals a substantial agreement of the observed luminosi- 
ties with model predictions. Indeed, observed bolometric 
magnitudes are mainly concentrated in the range from 
—4.5 to —6.4, while models find C stars for Mt, i < —4.7 
at solar metallicity, and for Mf, i < —4.5 at Z = Zq/3. 
Moreover, in commenting Figs. ^ a- n( i 01 we already no- 
ticed that the short wavelength filters are insufficient to 
estimate the magnitudes of red C-rich sources, enshrouded 
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by dust. C-star bolometr ic magnitudes below M \ in i = — 5 
do exist in LMC (see e.g. Ivan Loon et al"lll997l) . but they 
should not represent the average population. In view of 
the fact that optical C(N) sources always show C/O ra- 
tios close to 1, this can be easily understood in terms of 
evolution: subsequent phases give birth to redder objects, 
more C-rich (e.g. classes CV5— CV6 of BKR2), radiating 
their flux mainly in mid-IR and having a luminosity func- 
tion shifted at higher luminosities (see e.g. VL2). 

Compared to previous work on galactic C-rich giants 
by BKR1 & BKR2, our sample (with more homogeneous 
photometry and with a wider coverage of IR wavelengths) 
further emphasizes the importance of high luminosity C- 
stars. In the luminosity functions (LFs) of the above au- 
thors (see e.g. Figure 5 of BKR1, where the CV sources 
with the best distances were plotted) about 43% of CV3- 
CV6 sources are brighter than Mf, i = —5, while in our 
sample (Fig. |HJ these account for about 60% of the total 
sample; similarly 33% of CV3-CV6 sources in the quoted 
papers lay below Mf, i = —4, while in our data their frac- 
tional number is about 2%. On the whole, our distribution 
is displaced to higher luminosities due to a higher weight 
of observations at long wavelengths. 

Our sample of galactic stars integrates previ- 
ous IR studies on du st-enshrouded C stars in LMC 
ijvan Loon et all ll998t VLl.VL2l. These works were 
forcedly restricted to a relatively small number of sources: 
in this limit, their bright LF had little overlap with 
the fainter one of optically-selected C stars (VL1, VL2). 
However, other selection criteria lead to intermediate av- 
erage magnitudes (around M^oi = — 5: see lWhitelock et all 
2003|. Our Fig.|5]now shows that the LF of galactic C stars 
is in fact continuous, unique and quite wide. A comparison 
with the bottom panel shows that our global distribution 
looks pretty much as a superposition, with gaps filled, of 
LFs previously obtained for different samples of C stars 
in LMC. 

As mentioned previously, remaining uncertainties on 
bolometric magnitudes for the reddest stars, for which we 
have sometimes an incomplete coverage of the IR flux, 
would further populate the high-luminosity tail of the his- 
togram, making our conclusions even stronger. The few 
points we have at low luminosity (M^oi > —4) correspond 
to the HC classes of BKR2. We refer to those papers for a 
detailed list of the possible physical mechanisms explain- 
ing them. 

A few sources in Fig. \7\ exceed model predictions for 
the luminosity. Let us discuss this in some detail . We recall 
that a well-defined relation exists on the AGB (IP aczvnski 
lffjlji : llben fc Renzinilll983l IWagenhuber & Oroenewegenl 
1998) between the core mass (mn) and the luminosity 
(L). Only if wh reaches the Chandrasekhar limit can an 
AGB star attain the (slightly model-dependent) limit of 
—7.1™. Values higher than this require an extra source of 
energy, as in the case of massive AGB stars (M > 5 Mq), 
which develop very high temperatures at the base of the 
co nvective envelop e, allowing the onset of hbb. In models 
bv lKarakasI (|2003), AGB stars with hbb can in fact reach 



-10 



hbb limit _ 
m u -L limit 



-6 h 



-4 - 



rf » " 



.1. 



K-[12.5] 

Fig. 9. The HR diagram of observed C stars adopting as a 
temperature indicator the near-to-mid IR colour K~[12.5], 
where contributions from both the photosphere and dust 
are present. 
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Fig. 10. The HR diagram of observed C stars in mid-IR, 
adopting the mid-IR colour [8.8]-[21.3] as abscissa. 



Mb i ~ —7.75. Normal AGB stars (i.e. without hbb) must 
instead remain well below the mjj — L limit, as toh never 
reaches the Chandrasekhar limit, due to efficient mass loss. 

This is even more true for C stars: very bright 
Intermediate Mass Stars (M = 5-9 Mq) can hardly be- 
come carbon rich, even in the absence of hbb, as any 12 C 
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dredged-up from He-burning layers is diluted over sev- 
eral solar masses o f oxygen-rich mate r ial in the convective 
envelope (see e.g . iBiiaso et ;TlfHHHl lOallino et ablll998l 
lAbia et alJl200lll2002l) . 

Recently, C s tars more mass ive than 4-5 Mq have 
been obtained bv lKarakasI l|200^ . but only for relatively 
low metallicities (Z = 0.004 to Z = 0.008). Such short- 
living stars should however be extinct by now. Formation 
of very bright C stars is in principle possible in the evo- 
lution of high-mass AGB objects (say, 6 to 9 Mq), when 
hbb ceases at the end, and if mass loss suitably limits 
the envelope mass, so that the last fe w thermal pulses 
make the photosphere become C rich llFrost et al.l Tl998: 
iKarakas fc Lattanziol Eool iKarakasl l2003(l . Statistically, 
however, these massive C stars should be extremely rare, 
both for their low weight in the initial mass function and 
for the very short duration of the theoretically permitted 
phases. We can further notice that in Fig. El all the points 
representing known AGB stars with reliable distance es- 
timates (large dots) stay (with one exception) below the 
limit for moderately massive stars (M < 5M Q ); points 
beyond this level can therefore be simply explained as 
due to erroneous classifications (C-rich supergiants, e.g. of 
Wolf-Rayet type). Alternatively, they might be AGB stars 
whose luminosity has been overestimated due to poor dis- 
tances or to uncertainties induced by variability. 

HR-like diagrams can be constructed also using mid 
IR colours that include the effects of dust. An example 
is in Fig. it shows that there are no semiregulars or 
irregulars for K— [12.5] > 3 and no Miras for K— [12.5] < 
2; large IR excesses are associated to Mira variability or 
to post-AGB stars, or to unknown objects near the AGB 
end. Probably, the spread in magnitudes for any colour is 
also an outcome of a spread in stellar masses. 

An interesting piece of information is added by Fig. 1101 
where the HR diagram is constructed using a purely mid- 
IR colour as abscissa. Here we are looking only at dust, 
photospheres should not be seen: we don't expect much 
further reddening on the emission from dust, which re- 
radiates the stellar flux. Indeed, the area covered by known 
long period variables is now quite small, defining a narrow 
sequence for sources along the AGB, and a limited colour 
extension of some more evolved sources, with a maximum 
colour for Miras of [8.8] — [21.3] = 2. All known stars redder 
than this are post-AGB objects; some unknown sources 
observed by MSX are present in the reddest area, and we 
suggest that they be classified as stars at the termination 
of the AGB stage or above it, in evolution toward the 
ejection of a planetary nebula. Large values of the mid-IR 
colours are in fact expected when the relatively warm dust 
of AGB stars becomes cooler, forming the extended (and 
often detached) shells of post-AGB stars. 

6. Mass Loss 

Consider now Fig. ^2 (panel a) ; it shows the relation be- 
tween mass loss rates and bolometric magnitudes. In panel 
a), referring to irregulars, semiregulars and Miras, there 



is a general increase of mass loss rates with luminosity, 
which cannot however be reduced to a simple analytical 
formula (as examples, power-law relationships with three 
representative exponents have been added). The fact is 
that mass loss rates and luminosity should be correlated 
as, while becoming more luminous, AGB stars become 
also increasingly cool, their pulsation becomes stronger 
and both facts should power more intense stellar winds. 
The spread in panel a), partially confusing the relation, 
should then be ascribed to the presence of other variables, 
possibly the mass, and/or the fact that different mecha- 
nisms power stellar winds in different evolutionary stages. 
This evolutionary hyp othesis was suggested by VL2 and 
Ivan Loon et alJ ((2005) as the cause of a spread similar to 
the one we show; they i n particular indicated that a switch 
from single scattering <Jjuralll984() to multiple scattering 
l)Gail fc Sedlmavr|ll98fil) of photons on dust grains as the 
star evolves can explain mass loss rates in excess of the 
momentum of radiation L/c. Notice that the points rep- 
resenting Miras stay in the upper part of the diagram, 
showing a steeper dependence of mass loss rates on lumi- 
nosity. A sharp steepening of mass loss rates, up to a final 
superwind phase, is in fact necessary to let the whole en- 
velope be ejected before the AGB termination. Including 
post-AGB stars and unknown sources in the plot, as done 
in Fig. 1111 panel b), introduces further complexities. We 
believe we have here two problems. On one side, as al- 
ready mentioned, the reddest objects radiate part of their 
flux at wavelengths longer than 45 /xm, so that the lumi- 
nosity we estimate for them is only a lower limit. On the 
other side, these sources are leaving (or have just left) the 
AGB and should have huge circumstellar envelopes with 
possibly detached shells escaping our spatial resolution. 
Any correlation is then hampered by poor spatial and/or 
spectral coverage. 

A good monitor of the mass loss rate is instead avail- 
able from the IR colours. This is illustrated in Fig. ^1 In 
particular, when mid-IR colours are used, especially if in- 
cluding the 21.3 /xm (with the characteristic emission fea- 
ture of C-rich post-AGB stars) , then the variation of mass 
loss rates with the stellar classes (and with evolution, if 
the hypothesis that Miras are on average quite evolved 
holds) is clearly distinguishable (panel a) . All known long 
period variables remain at values of the [8.8] — [21.3] /im 
colour below 2, while post-AGB objects display much red- 
der colours, probably due to cold dust formation. Again, 
transition objects exist among the sources of unknown 
classification. A similar effect is seen from the K— [21.3] 
colour (panel b), but here the separation of post-AGB 
stars from the rest is less clear, and a rather continuous 
distribution is seen. In general, the growth of the IR excess 
linking Miras to post-AGBs seems to be rather smooth, 
with no obvious break (e.g. for the superwind phase that 
is expected to terminate the AGB evolution). 

Data like those discussed in this paper might be useful 
to infer which part of the mass return to the interstellar 
medium is due to carbon stars. This would require to in- 
tegrate the mass lost by AGB stars over an Initial Mass 
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Fig. 11. Top panel: the bolometric magnitudes for AGB 
C-rich stars of known variability type, as indicated, plot- 
ted as a function of mass loss rates. Bottom panel: same 
as before, but including unknown and post-AGB sources. 
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Fig. 12. Relations linking the mass loss rates of C-rich 
AGB and post-AGB sources to the IR colours up to 21 
/im. 



7. Conclusions 

In this paper we have presented the first issue of a se- 
ries of works on the IR properties of galactic AGB and 
post-AGB stars, based on archived data from the 2MASS, 
MSX (6C Point Source Catalogue) and ISO-SWS experi- 
ments. In particular, we analyzed a sample of more than 
230 C-rich AGB stars, for which estimates of distances, 
mass loss rates and colours from 1 to 21 (and sometimes 
45) ^m could be collected. We found that uncertainties 
still affect distances, except for a limited subset of objects 
with astrometric-quality data. Many sources were found 
to be extremely red, with SEDs extending to wavelengths 
unreachable by ground-based observations: in particular, 
we showed that, starting from the Mira stage, AGB car- 
bon stars begin to radiate most of their flux beyond 10 /zm. 
This excludes the possibility of deriving meaningful results 
on AGB luminosities by extrapolating fluxes obtained at 
optical or ncar-IR wavelengths (a practice that is unfor- 
tunately common even today). Using the best group of 
data, for which reliable distances and the whole SEDs up 
to 45 /im were available, we derived precise bolometric cor- 
rections, to be applied to the other sources, and thus we 
constructed the colour-magnitude diagrams in near and 
mid IR. The results suggest that AGB carbon stars of 
the Galaxy are objects of moderate mass (below about 5 
Mq), and reach bolometric magnitudes between —4.5 and 
—6.4. These findings confirm (at a higher metallicity and 
with larger statistics) previous suggestions on Magellanic 
Clouds by other authors. The luminosities derived are in 
good agreement with those predicted by stellar models 
using the Schwarzschild's criterion for convection and are 
at odd with the belief (still persisting in some theoreti- 
cal studies and population-synthesis works), according to 
which observed C-star luminosities would be lower than 
predicted by canonical stellar evolution. 

We also analyzed mass loss rates and their relation 
with the magnitudes: we suggested that correlations be- 
tween mass loss and luminosity are hard to derive for 
post-AGB objects, probably due to insufficient spatial and 
spectral coverage of their emission. For AGB stars in- 
stead, correlations of mass loss rates with luminosity ex- 
ist, but reveal dependence by more than one parameter. 
Here the stellar mass and the activation of different mech- 
anisms for powering mass loss should add their effects. 
Clearer correlations emerged between the mass loss rates 
and the IR colours, especially when the baseline was suf- 
ficiently extended in wavelength. These correlations also 
provided tools for classifying the evolutionary status of 
carbon stars. 



Function and a Star Formation Rate for low and interme- 
diate mass stars and to estimate the fractional contribu- 
tion of carbon-rich phases. We shall therefore reconsider 
this task in a subsequent issue of this analysis, after M 
giants are included in the discussion. 
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Table 1. Relevant data of our sample of carbon stars with astrometric, or reliable, distance estimates. Here and in 
next tables the following rules apply for labelling, i) For the variability type, M means Mira; S means Semiregular; 
I means Irregular (P is used for all Post-AGB star s, "-" for all stars of unknow n variabil ity), ii) References for mass 
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^kinnCT^taLn[ l ^^|) • ,Iaml991 for I.Iaming LeLall dlQQlf) . h i ) Refer ences f or distance are: B fo r | B^rgeatfc Ghevalherl 



L for Loud et al. 



Skinnci 

ESI 



G for iGroenewegen et alJ l|2002bh : S for lSchoier & Olofssonl J2OO1I). u sing methods 
fromlGroenewegen et all l)2fl02bfk S* for the same authors, w hen using methods fromlLoup et al.l 11199.^ : Meil998 for 
iMeixner et all l|l998|) ; Men2002 forFMen'shchikov et alJ l|2002l): S1997 fo r lSkinner et al.l <ll997ft: Hon200 3 for lHonv et all 
2003J); Bains2003 for lBains et all l)2003|) : Philll998 for lplmlips] l)l998|) : K: method bv lKastnerl (|l992|) . 



IRAS 


Other 


Var. 


Mass Loss 


Type 




d 


Type 


J 


H 


K 


8.8 


11.7 


12.5 


14.6 


21.3 


Data 


name 


name 


type 


[Mn/yrj 

L (1)/ J J 


ref. 


[km/s] 


[kpc] 


ref. 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


[Jy] 


origin 


09070 9fi 1 

uzz iu — zoiy 




M 


t. oni? nfi 


yj 


1 fi f. 

1D.0 


n sfi 




32 4 


QO O 

yy . u 


1 Q'i 
Do 


977 
Z ( 1 


209 


1 Qfi 
DO 


133 

loo 


SO -1 
oU. 1 


IO*. J 


UoZzy+4 Pzl 


V 0o4 jr CI 


1V1 


rr nnp nfi 


[5 

U 


15.0 


().i2 


"D 

rs 


21.0 


93. 1 


231 


lol 


52 1 


373 


204 


1!)!) 


1 Q 
1M ) 


ClA Kf\ I _l_ 10-10 
U4oU4+4y4y 


ATT A,,. 
AU A 111 


1V1 


c vnT? fl7 


u 




1.47 


[3 

L) 


13.2 


27. 5 


45.0 


16.4 




10. ( 


6. 1 


4.2 


MQY 
IVlOA 


U01o0-|-oZZ i 


UV Aur 


M 


1 finp nfi 




1 n n 


1 on 
i.uy 


JJ 


1Q O 


fi"3 K 
O.i.O 


yo.o 


",0 K 

ou.o 




37 fi 
■> 1 .0 


99 9 


1 -3 -3 
lo.o 


MSX 


0fi3 19J-039S 
U0o4Z-|-UoZo 


V Doo ivion 


J\I 


9 /ini? nc; 




l.i.l 


1 "?7 


JJ 


97 
U.z ; 


9 r; 
z.o 


19 () 


l OP; 
DO 




1 (Yi 
DO 


ftft ft 

nn.o 


OO 

ou.u 


MSX 


OQnftft qoi') 

U0U00 — oZ4J 


■\/0,1fi "P1111 

VJ40 rllp 


1V1 


Q 'i 1 T-P n 

O.OlJ-J-UO 


L> 


20.7 


1.12 


Li 


1.5 


8.2 


33.7 


177 




144 


141.2 


93.1 


T\,TQ Y 


Uy40Z+LooU 


w J_jeo 


1V1 


3.30E-05 


Li 


14.7 


0.15 


[3 

Li 


2 . i 


74 . 9 


409 


3o3 32 


3 / 100 


335 lo 


9fi'lOO 

zo.)Uy 


18191 


ren 
lo\J 


1 9 1 A 7-1-0/1 9^ 
1Z44 I -|-U4Z0 


T.TT \Fir 
r\XJ V 11 


J\I 


9 imf nfi 

Z. JUJ-j-UD 


jj 


1 fi 1 

In. 1 


O ('V>i 
U.Oo 


jj 


y. ( 


27 2 


82 


10 / 


1 79 

1 1 z 


1 "i 3 

lOO 


01 fi 
y 1.0 


^-i 3 

04. 


rcpi 
loi. J 


1 177_l_3Q/13 
104 1 I -|->>y40 


V PrR 

V LlD 


J\I 


1 *.nTT nfi 


jj 




n H/i 
1 1 . 1 


jj 


fi^ 

OO.U 


1 3 1 

10 1 


1 Oft 

1 jo 


on 
yu . u 


QQ fi 

j y.u 


fi9 cr 
OZ.O 


37 7 


93 ^ 
ZO.O 


1 ft; ri 

lOl. J 


1 1 OOO+Ool 


T Dra 


1\1 


9nT? nfi 




1 3 . 5 


1 .43 


[3 

L) 


37.8 


107 


1!)0 


220 


2 13 


15!) 


lOo 


07.0 


1 Q P 1 
lov I 


9n30f?_l_,1 7^7 
ZUoyO+4 (0 1 


V Cyg 


1V1 


fi '-inTP nfi 




11.5 


. 74 


13 

L) 


02. 1 


317 


.1!)!) 


i 08 


1 1 1 


512 


3 15 


22s 


1M ) 


91 390j_3fi^o 

Zl 6£\)-\-ooo\) 


V 14ZD oyg 


J\I 


/inp nfi 




1 1 
1 l.U 


n S9 




71 Q 

1 1 . y 


lyo 


q 1 n 


9 '"ift 




1 O - ' 

iyo 


IOU. ( 


1 03 3 


MSX 


91 3^H_l_*7K9^i 
ZL60o^r 1 oZJ 


tl Pirn 

^ep 


J\l 


9 Qni? nfi 




2'? [) 


<~,i\ 




1 7vj 
±10 


10^ 
lU.i 


p;ft", 

OoO 


lOO 
l:)U 


A ^ 
■loo 


• >ZO 


215 


121 


1 ft; n 

10 V J 


ZJoZU+4,510 


L/Ir^ AUG 


1V1 


9 1 ni-P n 


L) 


14.0 


0.84 


Li 


0.23 


2.9 


l!). 1 


o99 


i 50 


0:11 


510 


AAO 
-L-iZ 


IQPi 

101.J 


U1Z40 — Jz4o 


R Scl 


S 


<\ inu nt; 


B 


16.9 


0.47 


B 




540 


7 13 


172 


136 


93.5 


03.1 


50.9 


ISO 


03 3 7 '1 _L fi990 
UoJ f 4-|-DZZy 


TT wi 


c 



9 nnir nfi 




20.(1 


0.53 


[3 

L) 




385 


41 y 


120 


119 


85.6 


54.9 


:12.0 


lo( J 


no ■>"7"7 1 rr 1 on 
033 f I +012U 


V400 r^Cl 


s 


9 Ar\J? (17 


B 


9.0 


0.53 


B 


132 


2()0 


308 


09.3 




45.0 


28.4 


15.0 


MSX 


041zV + o030 


SY Pei 


s 


1 z.r\T? nfi 


B 


17.5 


1.43 


B 


23. 1 


62.7 


93.0 


30.8 




22.5 


10.8 


0.3 


MSX 


U44yT + ooZ0 


VJ40 AU1 


s 


n (mil 1 n"7 


B 


7„ 


1.11 


B 


91.7 


li)2 


185 


20.8 




14.3 


7.1 


5.0 


MSX 


ncLn9Q_Lm nfi 
UoUzo+UlUO 


V\ Ull 


Q 




q onT-P fl7 
O.SUJ-j-U . 


L> 


11.0 


0.41 


Li 


48!) 


941 


1028 


190 


180 


121 


80. 


43.1 


iftPi 
101.J 


ni^oisx'! a nfi 
U0Z^o+o4U0 


O Alll 


c 



/1 9nT? nfi 
4.ZU1-J-UD 


[3 


25.o 


1.13 


L> 


1 / . 8 


5 < . 1 


112 


123 




9:1.5 


58.4 


34.3 


MO A 


n^191 1 Ol 9/1 
U04Z1-|-Z4Z4 


rpTT rp 

1 u iau 


C 
O 


rr 7ni? f\7 


[3 




1.05 


[3 

L) 


74.1 


149 


150 


33. 


32. 1 


28.1 


19.7 


9.9 


lo( J 


nrr 1 1 ori/in 
054ZO+2U4U 


Y Tan 


s 


1 finTP nc 


B 


11.0 


0.74 


B 


242 


495 


181 


130 




77.8 


59.1 


36.4 


MSX 


0fin77_i_9fini 


TTT Pom 

1 u oein 


c 



1 IOI- n 07 
4.4U1L-U I 




11^ 
11.0 


n p;9 
u.oz 




1 7 1 
1 1 1 


'! 9 7 

oz / 


31 ^ 


7fi 7 
1 O. 1 




00.4 


90 O 


10.0 


MS X. 


06192+0722 


BN Mon 


g 


4.20E-06 


B 


24.2 


1.28 


13 


24.9 


59.1 


78.6 


30.4 




25.4 


11.2 


7.1 


MSX 


ofii 07_i_n397 
uoiy ( +uoz ( 


r u ivion 


Q 
O 


9 finf? c\7 


[3 

L> 


~ 


1 .45 


Li 


09.9 


148 


113 


16.2 




10. (i 


3.8 


3.7 


\ 1 Q Y 
A lo A 


nfi3i 1 fine 
Uoolo+lOUO 


Un. Lrem 


Q 
O 


q mTT 07 


[3 

L> 


14.9 


0.02 


[3 

Li 


o!).2 


158 


102 


AA fi 
44.0 




36.7 


15.6 


14.2 


Pi 1 Q V 
1\10A 


ofi^fi-LOfii -1 

U0000-|-U014 


iv.v ivion 


C 
O 


1 7np n7 

1. PUJ-j-U . 


JJ 




n fi7 


jj 


"iO 

oy.u 


11'-* 


l'-'ii 


9fi K 
ZO.O 




qi ^ 

ol.O 


1 1 7 
11. ; 


fi 9 
o.z 


MS K 


ncrroc ni 1 n 
UDOoO — UolU 


■\/"fi1 /I ATcin 

VD14 IvlOIl 


c 



1 Km? ns 


[3 




0.48 


[3 

L) 


8 1 . 6 


122 


131 


20.8 




10.5 


9.0 


3.(i 


i\10A 


n7H A K H79Q 
U (U40 — U ( Zo 


it i ivion 


c 



rr oriP f\7 


[3 


11.0 


O.tiO 


L> 


132 


2->^ 1 


2 ? j 


00.2 




44 . ( j 


21.0 


1 3 . 5 


IVIOA 


Oft ^oc;_l 1 79^ 
UoOZO+1 (ZD 


a unc 


c 



fi 9ni? c\7 


[3 


1 .0 


0.71 


[3 

L) 


3 7(> 


()■}.') 


()2!) 
















y nya 


c 



- . 7np n7 





y.u 


n !Q 




9 on 
zuu 


IRA 
.nil 


/Ii "i 


OO r i 




52 1 


3^ 1 
00.1 


I 7 

I I .y 


MSX 


00^89 ^Q^ft 

uyooz— oifoo 


tJ7 Par 


c 



1 finT-T n7 


JJ 


14 


O '"i7 




o.J.Z 


146 


1 ^3 
loo 


'"iO 1 
oU. 1 




9Q 


Q 9 

y.z 


p; 


MSX 


1 999fi_l.ni 0.9 
IZZZO-rUlUZ 


Oo V 11 


c 



*. RnT-P n7 


JJ 


£9 


O X.K 


jj 


81 Q 

oi.y 


L'.)-i 


■"i 1 n 
. > 1 U 


1 IV! 

DO 


1UD 


02 2 


3 1 1 
1. 4 


24 6 


i ft; pi 

lOV J 


19-1 07 _l -1 ^ 1 9 
1Z4Z ( +4D4Z 


1 UVll 


Q 
O 


1 mT-P n7 


[3 
L) 


8.5 


0.20 


[3 

Li 


(ill 


1331 


131 


2()'> 


102 


85.7 


53.8 


3.''i. 5 


rcn 


1 9^/1 1_i_fifi1 r. 
lZ044-r~0010 


RY Dra 


c 
i_> 


1 mi-P n7 


JJ 


1 n n 
111.11 


U.OO 


JJ 


"■*18 


425 


10 1 


1 97 
1Z 1 


yj.D 


73 1 
( 0.1 


0U.4 


9-1 7 
Z4. / 


1 ft; Pi 

lOl. J 


1 ^OQfi fiOOO 

iouyo— ouuy 


AO Vjll 


O 


7 ant? n7 






O ( Yi 

u.yo 


g 


",fi fi 
00.0 


1t> 1 


1 fi7 
10 1 


P. 9 r i 
OZ.O 




44 9 


21 


1 O O 

iu . y 


MSX 


1 71 -1 1 3^-11 
1 i 441 — 0041 


OA OCO 


c 



- . snp n7 






n s - ^ 

U.00 


g 


91 4 


1 a 1 

lol 


1 ft", 
loO 


37 fi 




25 4 


11 2 


7 1 
1.1 


MSX 


1 8476 — 0758 


S Set 


s 


5 60F-06 


13 


17.3 


0.58 


13 


191 


358 


371 


63.1 


53.5 


48.2 


28.0 


18.0 


ISO 


10.ni 7 n^il^. 

iyui 1 — uo4o 


V Aql 


c; 


fi finp n7 


13 


8 5 


56 


13 


■">s 1 


575 


71 


144 


117 


93 


66 1 


42 2 


ISO 


1Q3Q0-L399Q 

iyoyu+ozzy 


TT Over 




1 9fl17 n7 


13 


13 5 


0'-* 


13 


^9 


99 4 


110 


19 1 


12 8 


10 9 


7 5 


4 8 


ISO 


901 1 ^_L383-1 
ZU110-|-0o04 


n.o L.yg 


c 



9 onT-P n7 


JJ 




n fi^ 


jj 


83 


lOO 


1 Oft 
Do 


A 1 ft 
41. 




9 1 


14. 


8 9 

o.Z 


MSX 


911fift 4^1/1 
Zl 1UO — 1 ±OX c * 


T Ind 




1 7nT-P n7 


13 





05 


13 


'^57 


119 


397 


53 7 


38 7 


30 7 


22 9 


12 5 


ISO 


91 3QQ_I_3M fi 

zioyy^ooio 


VzlfiO rVn- 
V 4DU wyg 




a ^nP 1 n7 


13 


10 


64 


13 


303 


599 


520 


84 


67 9 


53 5 


31 


19 1 


ISO 


9".^S7J-finO_l 




c; 


1 - .nTP n7 


13 


3 8 


65 


13 


182 


330 


353 


53 5 


30 1 


23 4 


18 3 


13 7 


ISO 


03075+5742 


V623 Caa 




1.10E-07 


B 




0.51 


B 


114 


213 


219 


43.1 


29.3 


25.4 


19.4 


7.5 


ISO 


05449+3036 


FU Aur 




2.80E-07 


B 




1.29 


B 


49.5 


92.0 


96.4 


17.0 




12.8 


6.2 


3.1 


MSX 


06225+1445 


BL Ori 




2.50E-07 


B 


9.0 


0.60 


B 


212 


105 


328 


52.7 




32.1 


18.6 


10.0 


MSX 


07057-1150 


W CMa 




6.10E-07 


B 


10.5 


0.79 


B 


178 


298 


295 


53.6 




35.8 


23.5 


13.2 


MSX 


12465-6129 


RXCru 




6.20E-07 


B 




0.85 


B 


65.7 


149 


182 


52.1 




36.7 


18.7 


11.3 


MSX 


17172-4020 


V1079 Sco 




1.20E-06 


B 




0.92 


B 


88.2 


214 


20 1 


58.2 


86.7 


54.8 


30.1 


20.7 


ISO 


18040-0941 


FX Ser 




3.10E-05 


B 


28.4 


1.23 


B 


2.2 


13.3 


55.4 














18155-1327 


ES Ser 




3.00E-07 


B 




0.91 


B 


53.1 


119 


145 


40.2 




27.1 


15.0 


7.8 


MSX 


18306+3657 


T Lyr 




8.00E-07 


B 


11.5 


0.65 


B 


187 


130 


50 1 


121 


69.4 


41.9 


29.9 


21.0 


ISO 


18448+0523 


DR Ser 




1.00E-06 


B 


20.1 


1.29 


B 


46.7 


104 


123 


23.9 




15.1 


7.7 


4.7 


MSX 


19147+2149 


CG Vul 




3.30E-07 


B 




0.81 


B 


11.1 


87.2 


112 


38.8 




23.9 


11.1 


7.6 


MSX 


21553+5015 


LWCyg 




3.50E-07 


B 




1.09 


B 


36.0 


81.3 


110 


31.1 




19.2 


9.3 


5.4 


MSX 


22423+6127 


DG Cep 




2.30E-07 


B 




0.80 


B 


48.8 


102 


119 


20.1 




18.0 


!).() 


3.8 


MSX 


23438+0312 


TX Psc 




3.20E-07 


B 


10.7 


0.32 


B 


529 


1038 


1065 


160 


103 


82.8 


60.4 


35.5 


ISO 


04395+3601 


RAFGL 618 


P 


2.00E-04 


Meil998 


19.5 


1.70 


Mcil998 


< 0.01 


0.03 


0.20 


193 


400 


450 


651 


1260 


ISO 


06176-1036 


Red Rectangle 


P 


1.00E-04 


Men2002 


5.0 


0.71 


Men2002 


3.7 


9.0 


23.0 


323 


3S0 


369 


377 


472 


ISO 




RAFGL 2688 


P 


7.00E-05 


S1997 


19.5 


1.20 


S1997 


< 0.01 


0.19 


0.19 


170 


619 


749 


1253 


3496 


ISO 




NGC 7027 


P 


1.60E-04 


Jaml991 


18.6 


0.65 


Bains2003 


0.23 


1.0 


1.6 


133 


357 


332 


117 


1005 


ISO 


21282+5050 




P 








2.00 


Meil998 


0.04 


0.05 


0.10 


27.9 


00.7 


01.4 


64.1 


64.8 


ISO 


21306+4422 


1C 5117 


P 








10.23 


PM111998 


0.05 


0.05 


0.11 


0.7 


12.9 


11.4 


14.3 


33.5 


ISO 



Guandalini et al.: Infrared Luminosities of C Stars, Online Material p 3 



Table 2. Relevant data for AGB C stars of known variability type with non-astrometric distances. 
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Table 3. Relevant data for AGB C stars of unknown variability type and with non-astrometric distances. First part. 
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ref. 
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Table 4. Relevant data for AGB C stars of unknown variability type and with non-astrometric distances. Second 
part. 
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